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Abstract 

We establish the existence of a conservative weak solution to the Cauchy 
problem for the nonlinear variational wave equation uu — c(u)(c(u)u x ) x = 0, for 
initial data of finite energy. Here c(-) is any smooth function with uniformly 
positive bounded values. 
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1 Introduction 

We are interested in the Cauchy problem 

u tt - c(u)(c(u)u x ) x = 0, (1.1) 

with initial data 

u(0, x) = Uq(x) , u t (0, x) = Ui(x) . (1.2) 

Throughout the following, we assume that c:f h R+ is a smooth, bounded, uniformly 
positive function. Even for smooth initial data, it is well known that the solution can 
lose regularity in finite time ([12]). It is thus of interest to study whether the solution 
can be extended beyond the time when a singularity appears. This is indeed the main 
concern of the present paper. In ([5]) we considered the related equation 

u t + f(u) x = ^J*f"(u)u 2 x dx (1.3) 

and constructed a semigroup of solutions, depending continuously on the initial data. 
Here we establish similar results for the nonlinear wave equation (1.1). By introducing 
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new sets of dependent and independent variables, we show that the solution to the 
Cauchy problem can be obtained as the fixed point of a contractive transformation. 
Our main result can be stated as follows. 

Theorem 1. Let c: Em [ft -1 , k] be a smooth function, for some k > 1. Assume that 
the initial data uq in (1.2) is absolutely continuous, and that (uq) x 6 L , tii e L 2 . 
Then the Cauchy problem (1.1)- (1.2) admits a weak solution u = u(t,x), defined for 
all (t, i) 6 E x E. In the t-x plane, the function u is locally Holder continuous with 
exponent 1/2. This solution t \— > u(t,-) is continuously differ entiable as a map with 
values in Lf oc; for all 1 < p < 2. Moreover, it is Lipschitz continuous w.r.t. the L 2 
distance, i.e. 

|K*,.)-«(s,.)|| La <L\t-s\ (1.4) 
for all t,s el. The equation (1.1) is satisfied in integral sense, i.e. 

(f> t u t — (c(u)<f>) x c(u)u x dxdt = (1.5) 

for all test functions <f> e C\. Concerning the initial conditions, the first equality in 
(1.2) is satisfied pointwise, while the second holds in Lf oc for p e [1,2[. 

Our constructive procedure yields solutions which depend continuously on the initial 
data. Moreover, the "energy" 

£{t) = \ j [ M ?M + c >(^)K(^)] dx (1-6) 

remains uniformly bounded. More precisely, one has 

Theorem 2. A family of weak solutions to the Cauchy problem (1.1)- (1.2) can be 
constructed with the following additional properties. For every t Gl one has 

S{t) < S = [uKx) + c 2 (u (x))(u )l(x^ dx. (1.7) 

Moreover, let a sequence of initial conditions satisfy 

\\{Uq) x - («o)x|| L 3 » IK - M lllL2 -»• > 

and Mq — > -u uniformly on compact sets, as n — > oo. T/ien one /ias ^/ie convergence of 
the corresponding solutions u n — > -u, uniformly on bounded subsets of the t-x plane. 

It appears in (1.7) that the total energy of our solutions may decrease in time. Yet, 
we emphasize that our solutions are conservative, in the following sense. 

Theorem 3. There exists a continuous family {/i t ; t G M} of positive Radon measures 
on the real line with the following properties. 
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(i) At every time t, one has /i*(M) = So- 

(ii) For eacht, the absolutely continuous part of \x t has density j(uf + c 2 ul) w.r.t. the 
Lebesgue measure. 

(Hi) For almost every i G 1, the singular part of fi t is concentrated on the set where 
c'(u) = 0. 

In other words, the total energy represented by the measure /i is conserved in time. 
Occasionally some of this energy is concentrated on a set of measure zero. At the times 
t when this happens, fi T has a non-trivial singular part and S(r) < Sq. The condition 
(iii) puts some restrictions on the set of such times r. In particular, if c'(u) 7^ for all 
u, then this set has measure zero. 

The paper is organized as follows. In the next two subsections we briefly discuss the 
physical motivations for the equation and recall some known results on its solutions. In 
Section 2 we introduce a new set of independent and dependent variables, and derive 
some identities valid for smooth solutions. We formulate a set of equations in the new 
variables which is equivalent to (1.1). Remarkably, in the new variables all singularities 
disappear: Smooth initial data lead to globally smooth solutions. In Section 3 we 
use a contractive transformation in a Banach space with a suitable weighted norm 
to show that there is a unique solution to the set of equations in the new variables, 
depending continuously on the data uo, u\. Going back to the original variables u, t, x, 
in Section 4 we establish the Holder continuity of these solutions u = u(t,x), and 
show that the integral equation (1.5) is satisfied. Moreover, in Section 5, we study 
the conservativeness of the solutions, establish the energy inequality and the Lipschitz 
continuity of the map t 1— > u(t, ■). This already yields a proof of Theorem 2. In Section 
6 we study the continuity of the maps t 1— > u x (t, ■), t 1— > u t (t, •), completing the proof 
of Theorem 1. The proof of Theorem 3 is given in Section 7. 

1.1 Physical background of the equation 

Equation (1.1) has several physical origins. In the context of nematic liquid crystals, it 
comes as follows. The mean orientation of the long molecules in a nematic liquid crystal 
is described by a director field of unit vectors, n e § 2 , the unit sphere. Associated with 
the director field n, there is the well-known Oseen-Franck potential energy density W 
given by 

W (n, Vn) = a |n x (V x n)| 2 + (3(V • n) 2 + 7 (n • V x n) 2 . (1.8) 

The positive constants a, (3, and 7 are elastic constants of the liquid crystal. For the 
special case a — (5 — 7, the potential energy density reduces to 

If (n,Vn) =a|Vn| 2 , 

which is the potential energy density used in harmonic maps into the sphere S> 2 . There 
are many studies on the constrained elliptic system of equations for n derived through 
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variational principles from the potential (1.8), and on the parabolic flow associated 
with it, see [3, 9, 10, 16, 22, 36] and references therein. In the regime in which inertia 
effects dominate viscosity, however, the propagation of the orientation waves in the 
director field may then be modeled by the least action principle (Saxton [29]) 

— J {<9 t n • d t n - W (n, Vn)} dxdt = 0, n ■ n = 1. (1.9) 

In the special case a — /3 — 7, this variational principle (1.9) yields the equation for 
harmonic wave maps from (l + 3)-dimensional Minkowski space into the two sphere, see 
[8, 31, 32] for example. For planar deformations depending on a single space variable 
x, the director field has the special form 

n = cos u(x, t)e x + sin u(x, t)e y , 

where the dependent variable u G M 1 measures the angle of the director field to the 
x-direction, and e x and e y are the coordinate vectors in the x and y directions, respec- 
tively. In this case, the variational principle (1.9) reduces to (1.1) with the wave speed 
c given specifically by 

c 2 (u) = a cos 2 u + (3 sin 2 u. (1-10) 

The equation (1.1) has interesting connections with long waves on a dipole chain in 
the continuum limit ([13], Zorski and Infeld [45], and Grundland and Infeld [14]), and 
in classical field theories and general relativity ([13]). We refer the interested reader to 
the article [13] for these connections. 

This equation (1.1) compares interestingly with other well-known equations, e. g. 

d 2 u~d x [p(d x u)} = 0, (1.11) 

where p(-) is a given function, considered by Lax [25], Klainerman and Majda [24], and 
Liu [28]. Second related equation is 

d 2 t u - c 2 (u)Au = (1.12) 

considered by Lindblad [27] , who established the global existence of smooth solutions 
of (1.12) with smooth, small, and spherically symmetric initial data in M 3 , where the 
large-time decay of solutions in high space dimensions is crucial. The multi-dimensional 
generalization of equation (1.1), 

d 2 u - c(«)V • (c(u)Vu) = 0, (1.13) 

contains a lower order term proportional to cc'|Vm| 2 , which (1.12) lacks. This lower 
order term is responsible for the blow-up in the derivatives of u. Finally, we note that 
equation (1.1) also looks related to the perturbed wave equation 

d 2 t u- Am + /(m,Vm,VVm) = 0, (1.14) 
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where f(u, Vw, VVm) satisfies an appropriate convexity condition (for example, f = u p 
or / = a(d t u) 2 + b\Vu\ 2 ) or some nullity condition. Blow-up for (1.14) with a convexity 
condition has been studied extensively, see [2, 11, 15, 20, 21, 26, 30, 33, 34] and Strauss 
[35] for more reference. Global existence and uniqueness of solutions to (1.14) with 
a nullity condition depend on the nullity structure and large time decay of solutions 
of the linear wave equation in higher dimensions (see Klainerman and Machedon [23] 
and references therein). Therefore (1.1) with the dependence of c(u) on u and the 
possibility of sign changes in c'{u) is familiar yet truly different. 

Equation (1.1) has interesting asymptotic uni-directional wave equations. Hunter 
and Saxton ([17]) derived the asymptotic equations 

(ut + u n u x ) x = ]^nu n -\u x f (1.15) 

for (1.1) via weakly nonlinear geometric optics. We mention that the x-derivative of 
equation (1.15) appears in the high-frequency limit of the variational principle for the 
Camassa-Holm equation ([1, 6, 7]), which arises in the theory of shallow water waves. 
A construction of global solutions to the Camassa-Holm equations, based on a similar 
variable transformation as in the present paper, will appear in [4] 

1.2 Known results 

In [18], Hunter and Zheng established the global existence of weak solutions to (1.15) 
(n — 1) with initial data of bounded variations. It has also been shown that the 
dissipative solutions are limits of vanishing viscosity. Equation (1.15) (n = 1) is also 
shown to be completely integrable ([19]). In [3 7]- [44], Ping Zhang and Zheng study the 
global existence, uniqueness, and regularity of the weak solutions to (1.15) (n = 1,2) 
with L 2 initial data, and special cases of (1.1). The study of the asymptotic equation 
has been very beneficial for both the blow-up result [12] and the current global existence 
result for the wave equation (1.1). 

2 Variable Transformations 

We start by deriving some identities valid for smooth solutions. Consider the variables 

R = u t + c{u)u x , « 2 ^ 

S = u t - c(u)u x , 

so that 

R + S R — S . . 

ut = —z—, u x = ^ — • (2-2) 

l zc 

By (1.1), the variables R, S satisfy 

Rt — cR T = -j-f-R 2 — S 2 ), 

(2.3) 

S t + cS x = £ c (S 2 -R 2 ). 
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Multiplying the first equation in (2.3) by R and the second one by S, we obtain balance 
laws for R 2 and S 2 , namely 



(R 2 ) t -(cR 2 ) x = f c (R 2 S-RS 2 ), 
(S 2 ) t + (cS 2 ) x = -£(R 2 S-RS 2 ) 
As a consequence, the following quantities are conserved: 



(2.4) 



E=-(u 2 + c 2 u 2 x ) = ^^ } M = -u t u x = ^—. (2.5) 



Indeed we have 

E t + (c 2 M) x = 0, 
M t + E x = 0. 



(2.6) 



One can think of R 2 /4 as the energy density of backward moving waves, and S 2 /4 as 
the energy density of forward moving waves. 

We observe that, if R, S satisfy (2.3) and u satisfies (2.1b), then the quantity 

F = R-S- 2cu x (2.7) 

provides solutions to the linear homogeneous equation 

F t -cF x = ^(R + S + 2cu x )F. (2.8) 

In particular, if F = at time t — 0, the same holds for all t > 0. Similarly, if R, S 
satisfy (2.3) and u satisfies (2.1a), then the quantity 

G = R + S - 2u t 

provides solutions to the linear homogeneous equation 

G t + cG x = Ur + S-2cu x )G. 
2c 

In particular, if G = at time t — 0, the same holds for all t > 0. We thus have 

Proposition 1. Any smooth solution of (1.1) provides a solution to (2.1)-(2.3). Con- 
versely, any smooth solution of (2.1b) and (2.3) (or (2.1a) and (2.3)) which satisfies 
(2.2b) (or (2.2a) ) at time t = provides a solution to (1.1). 



The main difficulty in the analysis of (1.1) is the possible breakdown of regularity 
of solutions. Indeed, even for smooth initial data, the quantities u x ,Ut can blow up 
in finite time. This is clear from the equations (2.3), where the right hand side grows 
quadratically, see ([12]) for handling change of signs of d and interaction between R 
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and S. To deal with possibly unbounded values of R, S, it is convenient to introduce 
a new set of dependent variables: 



w = 2 arctan R , 



so that 



R = tan 



w 



z = 2 arctan S , 



S = tan 



Using (2.3), we obtain the equations 

2 



w t - cw x = - (R t - cR x ) 



1 + i? 2 



1 + 



(5 t + c5 a 



c' fi 2 - 
2c 1 + i? 2 ' 



c ; ^ 2 - i? 2 
2c 1 + S 2 ' 



(2.9) 



(2.10) 



(2.11) 



X = const. 



Y = const. 




Figure 1: Characteristic curves 

To reduce the equation to a semilinear one, it is convenient to perform a further 
change of independent variables (fig. 1). Consider the equations for the forward and 
backward characteristics: 



x + = c(u) 



x = —c[u) . 



(2.12) 



The characteristics passing through the point (t, x) will be denoted by 

s x + (s,t,x) , s x~ (s,t,x) , 

respectively. As coordinates (X, Y) of a point (t, x) we shall use the quantities 

"X~ (0,t,x) 



X 



f 

Jo 



(1 + R 2 (0,x))dx : 



I (1 + S 2 (0,x))dx. (2.13) 

'x+(0,t,x) 
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Of course this implies 

X t -c(u)X x = 0, Y t + c(u)Y x = 0, (2.14) 

(X x ) t - (cX x ) x = , (Y x ) t + (cY x ) x = . (2.15) 



Notice that 



■y rx (0,t,x+h) 

X x (t, x) = lim - / (1 + R 2 (0, x)) dx , 

h J x -(o,t,x) 

Y x (t, x) = lim - / (1 + ,S 2 (0, x))dx. 

h^O h J x +( ,t, x ) 



'x+(0,t,x) 

For any smooth function /, using (2.14) one finds 

ft + cf x = f x X t + f Y Y t + cf x X x + cj Y Yx = (X t + cX x )f x = 2cX x f x , 
ft ~ cf x = f x X t + f Y Y t - cf x X x - cf Y Y x = (Y t - cY x )f Y = -2cY x f Y . 

We now introduce the further variables 



X x Y x 



q ' 2c 1 + R 2 ' " p A 2c 1 + S 2 

Therefore 

& w w = i ( si * 2 f cos2 1 - sin2 ! cos2 f ) 9 > 

6 W TW = 6 ( sin2 f cos2 f - sin2 I cos2 f ) P • 




(2.16) 



P= v , <1 = — T^- (2.17) 



Notice that the above definitions imply 

(Xx)- 1 = ^ = P cos 2 I , {-Y x r = ^ = ? cos 2 Z - . (2.18) 

From (2.10)-(2.11), using (2.16)-(2.18), we obtain 

n (1 + S 2 ) c' i? 2 - S 2 (1 + R 2 ) d S 2 - i? 2 
2c^ -w Y = — , 2c^ >-z x 



(2.19) 



Using trigonometric formulas, the above expressions can be further simplified as 

w y = ^2 (cos z — cos w) q , 
zx = (cos w cos z) p . 
Concerning the quantities p, q, we observe that 

c x = c'u x = c'^p. (2.20) 



Using again (2.18) and (2.15) we compute 

Pt - c Px = (Xx)- 1 2R(R t - cR x ) - (X x )- 2 [(X x ) t - c(X x ) x ](l + R 2 ) 
= {X x y l 2R^ c {R 2 - S 2 ) - (X x )- 2 [c x X x ](l + R 2 ) 

= -s 2 )i- ^£(R -s)(i + r 2 ) 

= f cT f w [S(l + R 2 )-R(l + S 2 )]. 

q t + cq x = (-K,)- 1 2S(S t - cS x ) - (-Y x )- 2 [(-Y x ) t + c(-Y x ) x ](l + S 2 ) 

= (-Y x )-i 2S £(S 2 - R 2 ) - {-Y x )~ 2 [c x {Y x )]{\ + S 2 ) 

= I+V 2 ^ 2 - R 2 )i - ^(S -R)(l + S 2 ) 

= ij^[R(l + S 2 )-S(l + R 2 )}. 
In turn, this yields 

Py = (Pt ~ cp x ) 1 ^T ) = (p t - c Px )± 

_ c ' s(i+R 2 )-fi(i+s 2 ) _ c ' r s hi (221) 

~ Ac 2 (l+R 2 )(l+5 2 ) — 4c 2 ll+S 2 1+R 2 \ PV y^-^ 1 ) 

= £ [tan f cos 2 | - tan f cos 2 f ] pq = £ si^i^ pq , 

qx = {qt + cq x )^ 7) = {q t + cq x )j- c ^ 

_ d r(i+s 2 )-s(i+r 2 ) d r r s_] / 2 22) 

_ 4c 2 (1+5 2 )(1+ J R 2 ) PI - 4c 2 L+/J 2 1+S 2 J PI (Z.ZZj 

= ^ [tan f cos 2 f - tan § cos 2 f] pg = ^ ^f^ pg . 
Finally, by (2.16) we have 



ux = (u t + cu x ) ^j^z = j- cJ ^p= ± (tanf cos 2 f ) p, 
«y = K-^)^ = ^«=^ (tan f cos 2 f) g . 



(2.23) 



Starting with the nonlinear equation (1.1), using X, Y as independent variables we 
thus obtain a semilinear hyperbolic system with smooth coefficients for the variables 
u,w, z,p,q. Using some trigonometric identities, the set of equations (2.19), (2.21)- 
(2.22) and (2.23) can be rewritten as 

Wy = jfr (cos z — cos w) q , 

(2.24) 

zx = ^2 (cos w cos z) p , 
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Py = ^3 [ sin z - sin w] pq , 

(2.25) 

qx = ^2 [sinw - sinz] pq, 



ux = SJ ^P, 



4c 

(2.26) 

u Y = ^q- 

Remark 1. The function u can be determined by using either one of the equations in 
(2.26). One can easily check that the two equations are compatible, namely 

„, sinw J | cosw | sinm „, 

= g|^-3 [ — 2 sin w sin z + cos u> cos z — cos 2 u> — sin 2 w + sin w sin z] 

(2.27) 

= gf^ [cos(w - z) - l]p<? 
= Myx • 

Remark 2. We observe that the new system is invariant under translation by 2n in 
w and z. Actually, it would be more precise to work with the variables = e lw and 
z^ = e lz . However, for simplicity we shall use the variables w,z, keeping in mind that 
they range on the unit circle [— ir, ir] with endpoints identified. 

The system (2.24)-(2.26) must now be supplemented by non-characteristic boundary 
conditions, corresponding to (1.2). For this purpose, we observe that u Q ,U\ determine 
the initial values of the functions R, S at time t = 0. The line t = corresponds to a 
curve 7 in the (X, Y) plane, say 

Y = <p(X), iGl, 

where Y = <p(X) if and only if 

l>X PX 

X= (1 + R 2 (0,x))dx, Y = - (1 + S 2 (0,x))dx forsome igE. 
Jo Jo 

We can use the variable a; as a parameter along the curve 7. The assumptions u G H 1 , 
ui G L 2 imply R, S G L 2 ; to fix the ideas, let 

£o = ^J [R 2 (0, x) + S 2 (0, x)} dx<oo. (2.28) 

The two functions 

rx pO 

X(x)= (1 + R 2 (0,x))dx, Y(x)= (1 + S 2 (0,x))dx 

Jo Jx 
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are well defined and absolutely continuous. Clearly, X is strictly increasing while 
Y is strictly decreasing. Therefore, the map X \— > (p(X) is continuous and strictly 
decreasing. From (2.28) it follows 

\X + (p(X)\ < AS . (2.29) 

As (t,x) ranges over the domain [0,oo[ xR, the corresponding variables (X,Y) range 
over the set 

n + = {(X,Y); Y> V (X)}. (2.30) 

Along the curve 

1 = {{X,Y)- Y = tp(X)}cR 2 

parametrized by x i— > (X(x), Y(x)), we can thus assign the boundary data (w, z, p, q, u)E 
L°° defined by 

{t-:™(t1; fill: — (231 > 

We observe that the identity 

w z 

F = tan tan - - 2c(ii)u a . = (2.32) 

is identically satisfied along 7. A similar identity holds for G. 



3 Construction of integral solutions 

Aim of this section is to prove a global existence theorem for the system (2.24)-(2.26), 
describing the nonlinear wave equation in our transformed variables. 

Theorem 4. Let the assumptions in Theorem 1 hold. Then the corresponding problem 
(2.24)- (2.26) with boundary data (2.31) has a unique solution, defined for all (X,Y) G 
R 2 . 

In the following, we shall construct the solution on the domain Q + where Y > tp(X). 
On the complementary set Q~ where Y < tp(X), the solution can be constructed in an 
entirely similar way. 

Observing that all equations (2.24)-(2.26) have a locally Lipschitz continuous right 
hand side, the construction of a local solution as fixed point of a suitable integral 
transformation is straightforward. To make sure that this solution is actually defined 
on the whole domain Q + , one must establish a priori bounds, showing that p, q remain 
bounded on bounded sets. This is not immediately obvious from the equations (2.25), 
because the right hand sides have quadratic growth. 

The basic estimate can be derived as follows. Assume 

Co = sup 



c'(u) 



4c 2 (u) 



< 00 . 



(3.1) 
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From (2.25) it follows the identity 

qx + Py = . 

In turn, this implies that the differential form p dX — q dY has zero integral along every 
closed curve contained in In particular, for every (X, Y) G consider the closed 
curve £ (see fig. 2) consisting of: 

- the vertical segment joining (X, tp(X)) with (X,Y), 

- the horizontal segment joining {X,Y) with ((/? _1 (F), Y) 

- the portion of boundary 7 = {Y = ip(X)} joining ((/? _1 (F), Y) with (X, <p(X)). 



Figure 2: The closed curve E 
Integrating along S, recalling that p — q — 1 along 7 and then using (2.29), we obtain 




f* 100 p(X', Y) dX> + £ {x) q(X, Y>) dY' = X- <p~\Y) + Y - <p(X) 



(3.2) 



<2(\X\ + |y|+45 ). 
Using (3.1)-(3.2) in (2.25), since p, q > we obtain the a priori bounds 



p(X, Y) = exp {^ {x) ^ [sinz - sinw] q(X, Y') dY') 




(3.3) 



<exp{2C (|X| + |y|+4£ )}. 
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Similarly, 

q(X, Y) < exp {2C (|X| + \Y\ + 4£ )} . (3.4) 

Relying on (3.3)-(3.4), we now show that, on bounded sets in the X-Y plane, the 
solution of (2.24)-(2.26) with boundary conditions (2.31) can be obtained as the fixed 
point of a contractive transformation. For any given r > 0, consider the bounded 
domain 

a r = {(X,Y); Y>ip(X), X<r, Y < r} . 
Introduce the space of functions 

A r = {/ :Q r ^R; ||/||* = ess sup e^ x+Y) \f(X, Y)\ < oo} 
L (x,Y)en r J 

where k is a suitably large constant, to be determined later. For w,z,p,q,u G A r , 
consider the transformation T(w, z,p, q, u) = (w, z,p, q, it) defined by 

w(X, Y) = w(X, (p(X)) + ^gy (cos z - cosw) qdY, 

(3.5) 

z(X,Y) = z(<p-\Y),Y) + f* 1(Y) ^gy (cos^ - cos z)pdX, 
p(X, Y) = l + £ {x) 4^[smz- sin w] pq dY}, 
q(X, Y) = l + ^gy [ sin w - sin z] pq dX} , 

&(X,Y)=u(XMX))+ f Y ^QdY 

In (3.6), the quantities p, q are defined as 

p = mm{p, 2e 2Co d x l + l y l +4 ^}, g = min{g, 2 e 2C ^ x \ + \ Y ^} . (3.8) 

Notice that p = p, q = q as long as the a priori estimates (3.3)-(3.4) are satisfied. 
Moreover, if in the equations (2.24)- (2.26) the variables p, q are replaced with p, q, 
then the right hand sides become uniformly Lipschitz continuous on bounded sets in 
the X-Y plane. A straightforward computation now shows that the map T is a strict 
contraction on the space A r , provided that the constant k is chosen sufficiently big 
(depending on the function c and on r). 

Obviously, if r' > r, then the solution of (3.5)-(3.7) on Q r > also provides the solution 
to the same equations on Q r , when restricted to this smaller domain. Letting r — > oo, 
in the limit we thus obtain a unique solution (w, z,p, q, u) of (3.5)-(3.7), defined on the 
whole domain Q + . 

To prove that these functions satisfy the (2.24)-(2.26), we claim that p = p, q = q 
at every point (X, Y) e Vt + . The proof is by contradiction. If our claim does not hold, 
since the maps Y \— > p(X,Y), X \— > q(X,Y) are continuous, we can find some point 
(X*,Y*) E tt + such that 

p(X, Y) < 2e 2C ^ x \ + \ Y \ + ^ , q(X, Y) < 2e 2C ^ x \ + \ Y \ + ^ (3.9) 
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(3.6) 
(3.7) 



for all (A", y) eQ* = Q+n{X <X*, Y < Y*}, but either p(X*,Y*) > ^c (\x\+\y\+4£ ) 
or q(X*,Y*) > | e 2Co(|x|+|Y|+4f )_ By (3.9), we still have p = p, q = q restricted to SI*, 
hence the equations (2.24)-(2.26) and the a priori bounds (3.3)-(3.4) remain valid. In 
particular, these imply 

p(X*, y*) < e 2 Co(\X\ + \Y\+4£ ) ^ q ^ x *^ Y *} < e 2Co(|X|+|Y|+4£ ) ? 

reaching a contradiction. □ 

Remark 3. In the solution constructed above, the variables w, z may well grow outside 
the initial range ] — 7r, 7r[. This happens precisely when the quantities R, S become 
unbounded, i.e. when singularities arise. 

For future reference, we state a useful consequence of the above construction. 

Corollary 1. If the initial data Uq,Ui are smooth, then the solution (u,p,q,w, z) of 
(2.24)-(2.26), (2.31) is a smooth function of the variables (X,Y). Moreover, assume 
that a sequence of smooth functions (w^, u™) m >i satisfies 

uniformly on compact subsets o/R. Then one has the convergence of the corresponding 
solutions: 

(u m , p m , q m , w m , z m ) -> (u,p,q,w,z) 
uniformly on bounded subsets of the X-Y plane. 

We also remark that the equations (2.24)-(2.26) imply the conservation laws 

qx+PY = o, (f) x -(f) r = o. ("0) 

4 Weak solutions, in the original variables 

By expressing the solution u(X, Y) in terms of the original variables (t,x), we shall 
recover a solution of the Cauchy problem (1.1)-(1.2). This will provide a proof of 
Theorem 1. 

As a preliminary, we examine the regularity of the solution (w, w, z, p, q) constructed 
in the previous section. Since the initial data (uo) x and u\ are only assumed to be in L 2 , 
the functions w, z,p, q may well be discontinuous. More precisely, on bounded subsets 
of the X-Y plane, the equations (2.24)-(2.26) imply the following: 

- The functions w,p are Lipschitz continuous w.r.t. Y, measurable w.r.t. X. 

- The functions z, q are Lipschitz continuous w.r.t. X, measurable w.r.t. Y. 

- The function u is Lipschitz continuous w.r.t. both X and Y . 
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The map (X, Y) \— > (t, x) can be constructed as follows. Setting f — x, then / = t in 
the two equations at (2.16), we find 



c = 2cX x x x , ( 1 = 2cX x t x , 

-c = -2CY; x Y , [ 1 = -2cy; ty , 
respectively. Therefore, using (2.18) we obtain 

1 (l+cos uj) p 




2X X 4 

1 (l+cos z) q 

ZYx 4 

1 _ (l+cos tu)p 

2cXl ~~ 4c ' 

1 (l+cos z) q 

-2cY x ~ A~c • 



(4.1) 



(4.2) 



For future reference, we write here the partial derivatives of the inverse mapping, valid 
at points where w,z ^ —n. 



X. 
Y, 



2c 



x (l+cos w)p ' I { (l+cos w) p ' 

< (4.3) 

2 I y = 2c 

(l+cos z)q ' ^ * (l+cos z)q 



We can now recover the functions x = x(X, Y) by integrating one of the equations in 

(4.1) . Moreover, we can compute t = t(X,Y) by integrating one of the equations in 

(4.2) . A straightforward calculation shows that the two equations in (4.1) are equiv- 
alent: differentiating the first w.r.t. Y or the second w.r.t. X one obtains the same 
expression. 

(l+cos w) p Y psin wwy 



XXY = 



4 4 

■ ; - - lLl r sin z — sin w + sin(z — w)] = Xyx ■ 



~ 32c 2 

Similarly, the equivalence of the two equations in (4.2) is checked by 

txv-tyx = ( X -f)y+m x = lxxY-{^du Y + ^c'ux) 

= Yg^l [ sin z — sin w + sin(z — w)] 



g pq 

16 c 3 



[(1 + cosw) sinz — (1 + cosz) sinw] = 0. 



In order to define u as a function of the original variables t, x, we should formally 
invert the map (X,Y) \— > (t,x) and write u(t,x) = u(X(t,x), Y(t,x)). The fact 
that the above map may not be one-to-one does not cause any real difficulty Indeed, 
given (t*,x*), we can choose an arbitrary point (X*,Y*) such that t(X*,Y*) = t*, 
x(X*, Y*) = x*, and define u(t*,x*) = u(X*, Y*). To prove that the values of u do not 
depend on the choice of (X*,Y*), we proceed as follows. Assume that there are two 
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distinct points such that t(X 1 ,Y 1 ) = t(X 2 ,Y 2 ) = t* , x(X 1 ,Y 1 ) = x(X 2 ,Y 2 ) = x*. We 
consider two cases: 

Case 1: X x < X 2 , Y x < Y 2 . Consider the set 

F X . = [(X,Y); x(X,Y)<x*} 

and call dT x * its boundary. By (4.1), x is increasing with X and decreasing with 
Y. Hence, this boundary can be represented as the graph of a Lipschitz continuous 
function: X — Y = (f)(X + Y). We now construct the Lipschitz continuous curve 7 
(fig. 3a) consisting of 

- a horizontal segment joining (Xi,Yi) with a point A = (Xa,Ya) on dT x *, with 

Y A = Y U 

- a portion of the boundary dT x *, 

- a vertical segment joining (X 2 , Y 2 ) to a point B = (Xb, Yb) on dY x *, with Xb = X 2 . 

We can obtain a Lipschitz continuous parametrization of the curve 7 : [£1,^2] l— * ^ 2 
in terms of the parameter £ = X + Y . Observe that the map (X, Y) 1— > (t, x) is 
constant along 7. By (4.1)-(4.2) this implies (1 + cosw)X^ = (1 + cos z)Y^ = 0, hence 
sin it; • = sinz • = 0. We now compute 

u(X 2 , Y 2 )- u(X 1 , Y 1 ) = / 7 (u x rfX + My dy) 

proving our claim. 

Case 2: X x < X 2 , Y"i > Y 2 . In this case, we consider the set 

r> = {(X,Y); t(X,Y)<t*}, 

and construct a curve 7 connecting (X!,^) with (X 2 ,Y 2 ) as in fig. 3b. Details are 
entirely similar to Case 1. 

We now prove that the function u(t, x) = u(X(t, x), Y(t, x)) thus obtained is Holder 
continuous on bounded sets. Toward this goal, consider any characteristic curve, say 
t I— > x + (t), with x + = c(u). By construction, this is parametrized by the function 
X ^ (t(X,Y), x(X,Y)), for some fixed Y. Recalling (2.16), (2.14), (2.18) and (2.26), 
we compute 

£ [u t + c{u)u x f dt = J^(2cX x u x ) 2 (2X t )- 1 dX 

= Jj; 2c (p cos 2 f ^ (|;2sinfcosf) 2 rfX (4.4) 
<IxAJX < C T , 
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for some constant C T depending only on r. Similarly, integrating along any backward 
characteristics t i— > x~(t) we obtain 



f 

Jo 



[u t -c{u)u x f dt<C T . (4.5) 



Since the speed of characteristics is ±c(w), and c(u) is uniformly positive and bounded, 
the bounds (4.4)-(4.5) imply that the function u = u(t,x) is Holder continuous with 
exponent 1/2. In turn, this implies that all characteristic curves are C 1 with Holder 
continuous derivative. Still from (4.4)-(4.5) it follows that the functions R, S at (2.1) 
are square integrable on bounded subsets of the t-x plane. 

Finally, we prove that the function u provides a weak solution to the nonlinear wave 
equation (1.1). According to (1.5), we need to show that 

= // 4>t[( u t + cu x) + (u t - cu x )} - {c(u)4>) x [(u t + cu x ) - (u t - cu x )} dxdt 

= II [<f>t ~ ( c 0)d (u t + cu x ) dxdt + ff[(f> t + {c4>)x\ (ut ~ cu x ) dxdt (4.6) 

= // [(f) t - (c0) J R dxdt + H[(j)t + (c0) J S dxdt . 

By (2.16), this is equivalent to 

fj{-2cY^R + 2 CXAX S + «u x X. + * r Y.)* { S- S >}** = . (47) 

It will be convenient to express the double integral in (4.7) in terms of the variables 
X,Y. We notice that, by (2.18) and (2.14), 

pq 

dxdt= 2c(l + Ri)(l + S*) dXdY - 
Using (2.26) and the identities 



J- 2 w 1+cos w ( R sin w 

— CDS . , — 



1+R 2 wo 2 2' I 1+R 2 2 ' 

1 _ Pri o2 z_ 1+cos z j S _ sin z 

1+S 2 ~ L ' Ub 2 — 2 ' ^ 1+S 2 ~ 2 ' 
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(4.8) 



the double integral in (4.6) can thus be written as 

JJ {2c±^0 y R + 2c 1 -±f ( f> x S + c'(^p 1 -±f-^ q 1 -±f) cf>(S-R)} 



(4.9) 



P9 A YrlV 

2c(l+i? 2 )(l+5 2 ) ayval 

= Il{T^P^ + T^^x + c ^(m-m) ( f ) ( s - R )} dxdY 

(sin w sin z — sin w cos 2 | tan ^ — sin z cos 2 ^ tan §) j dXdY 
= //{ 0y + ^ 0x + ^ [ cos(w + z) - 1] 0} rfXrfF . 
Recalling (2.30), one finds 

(psjF) y + (£ s^ = (2c ^ )y + {2cUy)x 

= 4c'u x uy + 4:CUxy 

( 4 - 10 ) 

= ^ sinw sinz + |^ [cos(w - 2) - 1] 

= ^[cos(w + *)-l]. 

Together, (4.9) and (4.10) imply (4.7) and hence (4.6). This establishes the integral 
equation (1.5) for every test function <p e C\. 



5 Conserved quantities 



From the conservation laws (3.10) it follows that the 1-forms pdX — qdY and | dX + 
I dY are closed, hence their integrals along any closed curve in th X-Y plane vanish. 
From the conservation laws at (2.6), it follows that the 1-forms 

Edx - (c 2 M)dt, Mdx-Edt (5.1) 

are also closed. There is a simple correspondence. In fact 

Edx- (c 2 M) dt = V -dX -IdY --dx - M dx + E dt = dX + ^- dY - - dt. 
v ; 4 4 2 4c 4c 2 

Recalling (4.1)-(4.2), these can be written in terms of the X-Y coordinates as 

(±^plP dx _(l^ph dY , (5.2) 

(l-cosnQp (l-co S *) g 

8c 8c K J 
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respectively. Using (2.24)-(2.26), one easily checks that these forms are indeed closed: 



(1 — cos w) p 
8 



cpq 
64c2 



sin ^(1 — cosw) — sinto(l — cos z) 



(1 - cos z)q 
8 



1 — cos w) p 
8c 



cpq 
64c3 



sin(w + z) — (sinw + sin z) 



[1 — cos z) q 
8c 



x 
(5.4) 



In addition, we have the 1-forms 



^= (1 + cosw) ^x- (1 + cos ^y, 
4 4 

dt= {1 + cosw)p dx+ {1 + msz)q dY, 



4c 



4c 



(5.5) 
(5.6) 



which are obviously closed. 




Figure 4: 

The solutions u = u(X, Y) constructed in Section 3 are conservative, in the sense 
that the integral of the form (5.2) along every Lipschitz continuous, closed curve in the 
X-Y plane is zero. 

To prove the inequality (1.7), fix any r > 0. The case r < is identical. For a 
given r > arbitrarily large, define the set (fig. 4) 



r={(X,y); <t(X,Y) <r, X<r, Y < r} 

By construction, the the map (X, Y) 1— > (t, x) will act as follows: 

Ah(t,o), B^(r,b), C ^ (0, c) , D h-> (0, d) , 



(5.7) 
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(5.8) 



for some a < b and d < c. Integrating the 1-form (5.2) along the boundary of T we 
obtain 

Jab (1 ~ C °8 SW)P dX - (1 ~ c ° sz)(? dY 

= j dc (i^plz dx - Hz^fis rfy _ j da (i-c««,)p rfx _ ^ (i-caBz) g dy 

= £ I L 2 (0, x) + c 2 (w(0, x)) ^(0, x)j . 
On the other hand, using (5.5) we compute 

fa \ U t( T > X ) + c2 ( M ( T ' X )) U l( T > X )} dx 

-- f (l~cosii))p ,y r 



'ABn{cosu>^-l} 

< £ - 



ABn{cosz^-l} 



(i-co Sz)g dy (5.9) 



Notice that the last relation in (5.8) is satisfied as an equality, because at time t — 0, 
along the curve 70 the variables u>, z never assume the value — tt. Letting r — > +00 in 
(5.7), one has a — > —00, b — > +00. Therefore (5.8) and (5.9) together imply S(t) < So, 
proving (1.7). 




x+Kh 



Figure 5: Proving Lipschitz continuity 

We now prove the Lipschitz continuity of the map t 1— > u(t, •) in the L 2 distance. 
For this purpose, for any fixed time r, we let /i T = fj,~ + /i+ be the positive measure on 
the real line defined as follows. In the smooth case, 

fi-(]a,b[) = \ I R 2 (r,x)dx, fj,+(]a, b[) = \ f S 2 (r,x)dx. (5.10) 

To define ^ in the general case, let 7 T be the boundary of the set 

T T = {(X,Y); t(X,Y)<r}. (5.11) 
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Given any open interval ]a,b[, let A = (X A ,Y A ) and B = (X B ,Y B ) be the points on 
7 T such that 

x(A) = a , X P — Y P < X A — Y A for every point P G 7 T with x(P) < a , 
x{B) =b, Xp-Y P >X B -Y B for every point P G 7 T with x(P) > 6 . 



A i T (]o,6[)=^(]a,6[)+ A i+(]a,6[), (5.12) 



Then 
where 

MM[) = / (i^^AY = (i^i^. (5.13) 

Jab ° Jab 8 

Recalling the discussion at (5.1)-(5.2), it is clear that fi~ , ^ are bounded, positive 
measures, and /x T (M) = S , for all r. Moreover, by (5.10) and (2.5), 

pb pb -I 

/ c 2 u 2 x dx< / -(P 2 + ,S 2 )cfe <2/i(]a,6[). 

For any a < b, this yields the estimate 

/■ft 

|-u(r, b) — u(t, a)\ 2 < \b — a\ I u 2 (r,y)dy<2K 2 \b — a\fi T (]a,b[). (5-14) 

J a 

Next, for a given h > 0, y G R, we seek an estimate on the distance \u(r+h, y)—u(r, y)\. 
As in fig. 5, let "y T +h be the boundary of the set T T+ h, as in (5.11). Let P = (Px, Py) 
be the point on 7 r such that 

x(P) = y , Xpi — Ypi < Xp — Yp for every point P' G 7 r with x(P') < x , 

Similarly, let Q = (Qx, Qy) be the point on ^y T+ h such that 

x(Q) = V , Xq> - Y Q , < Xq-Yq for every point Q' G ^ T+h with x(Q') < y . 

Notice that X P < X Q and Y P < Y Q . Let P + = (X + ,Y + ) be a point on 7 T with 
X + = Xq, and let P~ = (X~,Y~) be a point on 7 r with F~ = Yq. Notice that 
x(P + ) E]y, y + , because the point (r, x(Q)) lies on some characteristic curve with 
speed — c(u) > —k, passing through the point (r+h, y). Similarly, x(P~) G ]y — nh, y[ . 
Recalling that the forms in (5.2) and (5.6) are closed, we obtain the estimate 

\u(Q) - u(P+)\ < \u Y (X Q , Y)\dY 



|y?|1N dY 



= I y^n± T ny ^=^y dY (515) 

< (jyt dY) V2 • ( j£ fi=Sl£ia d Y ) V2 
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The last term in (5.15) contains the integral of the 1-form at (5.2), along the curve 7 T , 
between P~ and P + . Recalling the definition (5.12)-(5.13) and the estimate (5.14), we 
obtain the bound 

\u(r + h,x) -u(t,x)\ 2 < 2\u(Q) - u(P + )\ 2 + 2\u(P + ) - u(P)\ 2 

< Ah • /i T (\x — Kh, x + Kh[) + 4/t 2 • (kH) ■ (J, T (]x , x + h[) . 

(5.16) 

Therefore, for any h > 0, 

\\u(t + h, •) - u(t, -)|| L 2 = (/ \u(t + h,x) — u(t, x)\ 2 dx) 1/2 

< (J 4(1 + K 3 )h ■ jj, T (}x - Kh, x + Kh[) dx) 1/2 

(5.17) 

= (4(/€ 3 + l)/ i VrW) 1/2 

= h-[4(^ + l)S ] 1/2 . 
This proves the uniform Lipschitz continuity of the map 1 1— > u(t, •), stated at (1.4). 



6 Regularity of trajectories 

In this section we prove the continuity of the functions t i— > u t (t, •) and t h- > u x (t, •), as 
functions with values in L p . This will complete the proof of Theorem I. 

We first consider the case where the initial data (uq) x , u\ are smooth with compact 
support. In this case, the solution u = u(X, Y) remains smooth on the entire X-Y 
plane. Fix a time r and let r y T be the boundary of the set T T , as in (5.11). We claim 
that 

l u (t,-)\ t=T = u t (r,-) (6.1) 
where, by (2.14), (2.18) and (2.26), 
u t (r, x) =u x X t + u Y Y t 

(6-2) 

p smw 2c |_ q smz 2c 



4c p(l+cos«>) 4c q(l+cosz) 2(l+cosui) 2(l+cosz) 

Notice that (6.2) defines the values of u t (r, ■) at almost every point x G K, i.e. at all 
points outside the support of the singular part of the measure \i T defined at (5.12). By 
the inequality (1.7), recalling that c{u) > re -1 , we obtain 

\u t (T,x)\ 2 dx < k 2 £(t) < k 2 £ . (6.3) 

Jr 

To prove (6.1), let any e > be given. There exist finitely many disjoint intervals 
[a*, bj\ C M., % — 1, . . . , N, with the following property. Call Ai,P>i the points on 7 T 
such that x(Ai) = a iy x(Bi) = 6j. Then one has 

min{l + cos w(P) , 1 + cos^(P)} < 2e (6.4) 
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at every point P on 7 T contained in one of the arcs AiB^ while 
1 + cos w(P) > e , 1 + cos z(P) > e , 



(6.5) 



for every point P along r j T1 not contained in any of the arcs A^Bi. Call J = Ui<j<Ar[cti, 
J' = R\J, and notice that, as a function of the original variables, u = u(t, x) is smooth 
in a neighborhood of the set {r} x J'. Using Minkowski's inequality and the differen- 
tiability of u on J', we can write 



lim^o \ (/ M w(t + /i, x) - u(r, x) - /i it t (r, x) rfxj 

p \ 1 /p 



i/p 



< lim 



u(r + h, x) — m(t, xj 



dxj + (J 7 |wt(r, x)| p dx) 



i/p 



(6.6) 



We now provide an estimate on the measure of the "bad" set J: 

meas (J) = Jjdx = il±^ rfX - fi±2»ia rfF 



<2e f 

— J-y. 



(1— cos u>)p 



(6.7) 



rfx _ (l-cos, )g rfy < 2g g . 



Now choose q = 2/(2 — p) so that | + ^ = 1. Using Holder's inequality with 



conjugate exponents 2/p and g, and recalling (5.17), we obtain 

p / 2 \ 

u(r + h, x) — u(t, x) dx < meas (J) 1 ^ ■ I Jj u(r + /i, x) — u(r, x) rfx J 

<[2sS ] 1/c! -[\\u(r + h,.)-u(r,.)\\l 2 ) 

<[2e£,] 1/q -(h 2 [A{^ + l)£ ]) 



\ P/2 



p/2 



Therefore, 



lim sup — 



u(t + h,x) — u(t, x) — h 
In a similar way we estimate 



dx 



i/p 



< [2^ ] 1 ^.[4( K 3 + 1)^ ] 1/2 . (6.8) 



2 \ P/2 

dx 



^ Mr,x)| p dx^ ^ < meas(J) 1/w • [k 2 £ ] p/2 . 



(6.9) 
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Since e > is arbitrary, from (6.6), (6.8) and (6.9) we conclude 



lim — 

h^O h 



u(t + h, x) — u(r, x) — hu t (r, x) 



dx 



i/p 



= 0. 



(6.10) 



The proof of continuity of the map t h- > u t is similar. Fix e > 0. Consider the intervals 
[dj, 6j] as before. Since u is smooth on a neighborhood of {t} x J', it suffices to estimate 

lim sup fe ^ f \u t (r + h, x) — u t (r, x) \ p dx 

< limsup^o J j \u t (r + h, x) — u t (r, x)\ p dx 



< limsup^o [meas (J)] 1 ^ • ^ Jj u t {r + h, x) — ti*(r, x) 

< limsup^o [2eS } 1/q ■ (\\u t (r + h, -)|| L 2 + |K(r, -)IIe^) 



\ P/2 



<[2^ ] 1/9 [4^of. 
Since e > is arbitrary, this proves continuity. 

To extend the result to general initial data, such that (u ) x ,Ui G L 2 , we consider 
a sequence of smooth initial data, with (uq) x ,u\ e C^°, with -Uq — >• uq uniformly, 



u, 



(uq) x almost everywhere and in L , u\ — > u\ almost everywhere and in L . 



The continuity of the function t h- > u x (t, •) as a map with values in L p , 1 < p < 2, 
is proved in an entirely similar way. 



7 Energy conservation 

This section is devoted to the proof of Theorem 3, stating that, in some sense, the total 
energy of the solution remains constant in time. 

A key tool in our analysis is the wave interaction potential, defined as 

A(t) = (pi; <g> /i+){(x, y) ; x > y} . (7.1) 

We recall that /if are the positive measures defined at (5.13). Notice that, if ^tiHt 
are absolutely continuous w.r.t. Lebesgue measure, so that (5.10) holds, then (7.1) is 
equivalent to 

m = \ [ [ R 2 (t,x)S 2 (t,y)dxdy. 

4 J Jx>y 

Lemma 1. The map t t— > A(t) has locally bounded variation. Indeed, there exists a 
one-sided Lipschitz constant L Q such that 

A(t)-A(s) <L -(t-s) t>s>0. (7.2) 
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To prove the lemma, we first give a formal argument, valid when the solution 
u = u(t,x) remains smooth. We first notice that (2.4) implies 

j t (4A(f)) <- J 2cR 2 S 2 dx + J (R 2 + S 2 ) dx ■ J ^-\R 2 S - RS 2 \ dx 



< -2k 



R 2 S 2 dx + 4£ 



c 



L°° 



\R 2 S - RS 2 \dx, 



where k 1 is a lower bound for c(u). For each e > we have \R\ < e x / 2 + e 1 / 2 R 2 . 
Choosing e > such that 

c 
2^ 



« _1 > 4£ n 



we thus obtain 



(4A(0) < 



-7 



16 £ 2 



c 

2^ 



L°° 



This yields the L 1 estimate 

jT J {\R 2 S\ + \RS 2 \) dxdt = 0(1) • [A(0) + S 2 r] = (9(1) • (1 + r)£ 2 , 

where 0(1) denotes a quantity whose absolute value admits a uniform bound, depend- 
ing only on the function c = c(u) and not on the particular solution under consideration. 
In particular, the map t i— > A(t) has bounded variation on any bounded interval. It 
can be discontinuous, with downward jumps. 

To achieve a rigorous proof of Lemma 1, we need to reproduce the above argument 
in terms of the variables X,Y. As a preliminary, we observe that for every e > there 
exists a constant k p such that 



|sinz(l — cosw) — siniu(l — cosz)\ 
< k £ ■ (tan 2 ^ + tan 2 |) (1 + cosiu)(l + cos z) + e(l — cosiu)(l — cos z) 



(7.3) 



for every pair of angles w, z. 

Now fix < s < t. Consider the sets T s ,T t as in (5.11) and define T st = T t \ T s . 
Observing that 



VQ 

dxdt = + cosw)(l + cosz)dXdY , 



we can now write 



III. 



s J—oo 4 



dxdt 



(t - s)£ 



JJ Tat \ (tan 2 f + tan 2 f) • f c (1 + cosw)(l + cos^) dXdY . 



(7.4) 
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The first identity holds only for smooth solutions, but the second one is always valid. 
Recalling (5.4) and (5.13), and then using (7.3)-(7.4), we obtain 

A(f) - A(s) < - ff Tst i=^p ■ l=f±qdXdY 

+S • JJ T f^zPQ [sinz(l — cosw) — sinu>(l — cos z)\ dXdY 

< — ^ JJ T (1 — cosu>)(l — cos z) pq dXdY 

+£ ■ JJ T ^zPq [« e • (tan 2 | + tan 2 f) (1 + cosiu)(l + cos z) 
+e(l — cosiu)(l — cos z)\ dXdY 

< n(t — s) , 

for a suitable constant k. This proves the lemma. 
To prove Theorem 3, consider the three sets 

n 1 = {(X,Y); w(X,Y) = -n, z(X,Y) ^ -n , ^(1,7))^}, 

n 2 = {(X,Y); z(X,Y) = -n, w(X,Y)^-n, c'(u(X, Y)) ^ j , 

{1 3 = {(X,Y); z(X,Y) = -n, w(X,Y) = -tt, c'(u(X, Y)) + o} . 

From the equations (2.24), it follows that 

meas (Sli) = meas (f^) = . (7.5) 

Indeed, wy ^ on Qi and ^ on 2 . 

Let be the set of Lebesgue points of Q 3 . We now show that 

meas ({t(X, Y) ; (X,Y) e = . (7.6) 

To prove (7.4), fix any P* = (X*, Y*) G fl* 3 and let r = t(P*). We claim that 

A(r - h) - A(r + k) 
hmsup — 1 = +oo (7.7) 

h,k^o+ h + k 

By assumption, for any e > arbitrarily small we can find 5 > with the following 
property. For any square Q centered at P* with side of length £ < 5, there exists a 
vertical segment a and a horizontal segment a', as in fig. 6, such that 

meas (fi 3 n a) > (1 - e)£ , meas (Q 3 n </) > (1 - e)£ , (7.8) 

Call 



t + = max Y) ; (X, F) G a U </} , 
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Figure 6: 

r = min Y) ; GdUdj. 

Notice that, by (4.2), 

t+-r< [ (1 + C A ° SW)P dX+ [ {1 + C ; Sz)q dY<c -(ei)\ (7.9) 

Indeed, the integrand functions are Lipschitz continuous. Moreover, they vanish oustide 
a set of measure el. On the other hand, 

A(r) - A(t+) > ci(i - £ ) 2 ^ 2 - c 2 (t + - r) (7.io) 

for some constant C\ > 0. Since e > was arbitrary, this implies (7.5). 

Recalling that the map 1 1 — > A has bounded variation, from (7.5) it follows (7.4). 

We now observe that the singular part of /i T is nontrivial only if the set 

{P G 7 T ; w(P) = -7r or z(P) = -vr} 

has positive 1-dimensional measure. By the previous analysis, restricted to the region 
where d ^ 0, this can happen only for a set of times having zero measure. 
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